Genomic and epigenomic analysis of plasma cell-free DNA identifies stemness features associated with worse survival in AR -altered lethal prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) resistant to androgen receptor (AR)-targeted agents is often lethal. Unfortunately, biomarkers for this deadly disease remain under investigation, and underpinning mechanisms are ill-understood. Here, we applied deep sequencing to ∼100 mCRPC patients prior to the initiation of first-line AR-targeted therapy, which detected AR /enhancer alterations in over a third of patients, which correlated with lethality. To delve into the mechanism underlying why these patients with cell-free AR /enhancer alterations developed more lethal prostate cancer, we next performed genome-wide cell-free DNA epigenomics. Strikingly, we found that binding sites for transcription factors associated with developmental stemness were nucleosomally more accessible. These results were corroborated using cell-free DNA methylation data, as well as tumor RNA sequencing from a held-out cohort of mCRPC patients. Thus, we validated the importance of AR /enhancer alterations as a prognostic biomarker in lethal mCRPC, and showed that the underlying mechanism for lethality involves reprogramming developmental states toward increased stemness.

Association between the location of social medical insurance and social integration among China's elderly rural migrants: a nationwide cross-sectional study.

BACKGROUND: Universal social medical insurance coverage is viewed as a major factor in promoting social integration, but insufficient evidence exists on the integration of elderly rural migrants (ERM), generally aged 60 years and above, in low- and middle-income countries. To address this problem, we explore the relationship between the location of social medical insurance (SMI), such as a host city, and social integration in the context of Chinese ERM. METHODS: This study is based on data from the 2017 National Internal Migrant Dynamic Monitoring Survey in China. The study participants were Chinese ERM. An integration index was constructed to measure the degree of social integration in a multi-dimensional manner using a factor analysis method. This study used descriptive statistics and one-way analysis of variance to explore the differences in social integration between ERM with SMI from host cities and hometowns. Stepwise multiple linear regression analysis was used to test the correlation between SMI location and social integration level in the overall sample. Finally, the results were verified by propensity score matching. RESULTS: It was found that 606 (18.2%) of the insured ERM chose host city SMI, while 2727 (81.8%) chose hometown SMI. The level of social integration was lower among ERM with hometown SMI (-1.438 ± 32.795, F = 28.311, p ≤ 0.01) than those with host city SMI (6.649 ± 34.383). Among the dimensions of social integration, social participation contributed more than other factors, with a contribution rate of 45.42%. Host city SMI increased the probability of the social integration index by 647% among ERM (k-nearest neighbor caliper matched (n = 4, caliper = 0.02), with a full sample ATT value of 6.47 (T = 5.32, SE = 1.48, p < 0.05)). CONCLUSIONS: ERM with host city SMI have a higher social integration level than those with hometowns SMI. That is, host city SMI positively affects social integration. Policymakers should focus on the access of host city SMI for ERM. Removing the threshold of host city SMI coverage for ERM can promote social integration.

PACT: a pipeline for analysis of circulating tumor DNA.

MOTIVATION: Detection of genomic alterations in circulating tumor DNA (ctDNA) is currently used for active clinical monitoring of cancer progression and treatment response. While methods for analysis of small mutations are more developed, strategies for detecting structural variants (SVs) in ctDNA are limited. Additionally, reproducibly calling small-scale mutations, copy number alterations, and SVs in ctDNA is challenging due to the lack to unified tools for these different classes of variants. RESULTS: We developed a unified pipeline for the analysis of ctDNA [Pipeline for the Analysis of ctDNA (PACT)] that accurately detects SVs and consistently outperformed similar tools when applied to simulated, cell line, and clinical data. We provide PACT in the form of a Common Workflow Language pipeline which can be run by popular workflow management systems in high-performance computing environments. AVAILABILITY AND IMPLEMENTATION: PACT is freely available at https://github.com/ChrisMaherLab/PACT.

Urine cell-free DNA multi-omics to detect MRD and predict survival in bladder cancer patients.

Circulating tumor DNA (ctDNA) sensitivity remains subpar for molecular residual disease (MRD) detection in bladder cancer patients. To remedy this problem, we focused on the biofluid most proximal to the disease, urine, and analyzed urine tumor DNA in 74 localized bladder cancer patients. We integrated ultra-low-pass whole genome sequencing (ULP-WGS) with urine cancer personalized profiling by deep sequencing (uCAPP-Seq) to achieve sensitive MRD detection and predict overall survival. Variant allele frequency, inferred tumor mutational burden, and copy number-derived tumor fraction levels in urine cell-free DNA (cfDNA) significantly predicted pathologic complete response status, far better than plasma ctDNA was able to. A random forest model incorporating these urine cfDNA-derived factors with leave-one-out cross-validation was 87% sensitive for predicting residual disease in reference to gold-standard surgical pathology. Both progression-free survival (HR = 3.00, p = 0.01) and overall survival (HR = 4.81, p = 0.009) were dramatically worse by Kaplan-Meier analysis for patients predicted by the model to have MRD, which was corroborated by Cox regression analysis. Additional survival analyses performed on muscle-invasive, neoadjuvant chemotherapy, and held-out validation subgroups corroborated these findings. In summary, we profiled urine samples from 74 patients with localized bladder cancer and used urine cfDNA multi-omics to detect MRD sensitively and predict survival accurately.

Correction: Urine tumor DNA detection of minimal residual disease in muscle-invasive bladder cancer treated with curative-intent radical cystectomy: A cohort study.

[This corrects the article DOI: 10.1371/journal.pmed.1003732.].

Cell-free DNA ultra-low-pass whole genome sequencing to distinguish malignant peripheral nerve sheath tumor (MPNST) from its benign precursor lesion: A cross-sectional study.

BACKGROUND: The leading cause of mortality for patients with the neurofibromatosis type 1 (NF1) cancer predisposition syndrome is the development of malignant peripheral nerve sheath tumor (MPNST), an aggressive soft tissue sarcoma. In the setting of NF1, this cancer type frequently arises from within its common and benign precursor, plexiform neurofibroma (PN). Transformation from PN to MPNST is challenging to diagnose due to difficulties in distinguishing cross-sectional imaging results and intralesional heterogeneity resulting in biopsy sampling errors. METHODS AND FINDINGS: This multi-institutional study from the National Cancer Institute and Washington University in St. Louis used fragment size analysis and ultra-low-pass whole genome sequencing (ULP-WGS) of plasma cell-free DNA (cfDNA) to distinguish between MPNST and PN in patients with NF1. Following in silico enrichment for short cfDNA fragments and copy number analysis to estimate the fraction of plasma cfDNA originating from tumor (tumor fraction), we developed a noninvasive classifier that differentiates MPNST from PN with 86% pretreatment accuracy (91% specificity, 75% sensitivity) and 89% accuracy on serial analysis (91% specificity, 83% sensitivity). Healthy controls without NF1 (participants = 16, plasma samples = 16), PN (participants = 23, plasma samples = 23), and MPNST (participants = 14, plasma samples = 46) cohorts showed significant differences in tumor fraction in plasma (P = 0.001) as well as cfDNA fragment length (P < 0.001) with MPNST samples harboring shorter fragments and being enriched for tumor-derived cfDNA relative to PN and healthy controls. No other covariates were significant on multivariate logistic regression. Mutational analysis demonstrated focal NF1 copy number loss in PN and MPNST patient plasma but not in healthy controls. Greater genomic instability including alterations associated with malignant transformation (focal copy number gains in chromosome arms 1q, 7p, 8q, 9q, and 17q; focal copy number losses in SUZ12, SMARCA2, CDKN2A/B, and chromosome arms 6p and 9p) was more prominently observed in MPNST plasma. Furthermore, the sum of longest tumor diameters (SLD) visualized by cross-sectional imaging correlated significantly with paired tumor fractions in plasma from MPNST patients (r = 0.39, P = 0.024). On serial analysis, tumor fraction levels in plasma dynamically correlated with treatment response to therapy and minimal residual disease (MRD) detection before relapse. Study limitations include a modest MPNST sample size despite accrual from 2 major referral centers for this rare malignancy, and lack of uniform treatment and imaging protocols representing a real-world cohort. CONCLUSIONS: Tumor fraction levels derived from cfDNA fragment size and copy number alteration analysis of plasma cfDNA using ULP-WGS significantly correlated with MPNST tumor burden, accurately distinguished MPNST from its benign PN precursor, and dynamically correlated with treatment response. In the future, our findings could form the basis for improved early cancer detection and monitoring in high-risk cancer-predisposed populations.

Urine tumor DNA detection of minimal residual disease in muscle-invasive bladder cancer treated with curative-intent radical cystectomy: A cohort study.

BACKGROUND: The standard of care treatment for muscle-invasive bladder cancer (MIBC) is radical cystectomy, which is typically preceded by neoadjuvant chemotherapy. However, the inability to assess minimal residual disease (MRD) noninvasively limits our ability to offer bladder-sparing treatment. Here, we sought to develop a liquid biopsy solution via urine tumor DNA (utDNA) analysis. METHODS AND FINDINGS: We applied urine Cancer Personalized Profiling by Deep Sequencing (uCAPP-Seq), a targeted next-generation sequencing (NGS) method for detecting utDNA, to urine cell-free DNA (cfDNA) samples acquired between April 2019 and November 2020 on the day of curative-intent radical cystectomy from 42 patients with localized bladder cancer. The average age of patients was 69 years (range: 50 to 86), of whom 76% (32/42) were male, 64% (27/42) were smokers, and 76% (32/42) had a confirmed diagnosis of MIBC. Among MIBC patients, 59% (19/32) received neoadjuvant chemotherapy. utDNA variant calling was performed noninvasively without prior sequencing of tumor tissue. The overall utDNA level for each patient was represented by the non-silent mutation with the highest variant allele fraction after removing germline variants. Urine was similarly analyzed from 15 healthy adults. utDNA analysis revealed a median utDNA level of 0% in healthy adults and 2.4% in bladder cancer patients. When patients were classified as those who had residual disease detected in their surgical sample (n = 16) compared to those who achieved a pathologic complete response (pCR; n = 26), median utDNA levels were 4.3% vs. 0%, respectively (p = 0.002). Using an optimal utDNA threshold to define MRD detection, positive utDNA MRD detection was highly correlated with the absence of pCR (p < 0.001) with a sensitivity of 81% and specificity of 81%. Leave-one-out cross-validation applied to the prediction of pathologic response based on utDNA MRD detection in our cohort yielded a highly significant accuracy of 81% (p = 0.007). Moreover, utDNA MRD-positive patients exhibited significantly worse progression-free survival (PFS; HR = 7.4; 95% CI: 1.4-38.9; p = 0.02) compared to utDNA MRD-negative patients. Concordance between urine- and tumor-derived mutations, determined in 5 MIBC patients, was 85%. Tumor mutational burden (TMB) in utDNA MRD-positive patients was inferred from the number of non-silent mutations detected in urine cfDNA by applying a linear relationship derived from The Cancer Genome Atlas (TCGA) whole exome sequencing of 409 MIBC tumors. We suggest that about 58% of these patients with high inferred TMB might have been candidates for treatment with early immune checkpoint blockade. Study limitations included an analysis restricted only to single-nucleotide variants (SNVs), survival differences diminished by surgery, and a low number of DNA damage response (DRR) mutations detected after neoadjuvant chemotherapy at the MRD time point. CONCLUSIONS: utDNA MRD detection prior to curative-intent radical cystectomy for bladder cancer correlated significantly with pathologic response, which may help select patients for bladder-sparing treatment. utDNA MRD detection also correlated significantly with PFS. Furthermore, utDNA can be used to noninvasively infer TMB, which could facilitate personalized immunotherapy for bladder cancer in the future.

ctDNA MRD Detection and Personalized Oncogenomic Analysis in Oligometastatic Colorectal Cancer From Plasma and Urine.

We hypothesized that circulating tumor DNA (ctDNA) molecular residual disease (MRD) analysis without prior mutational knowledge could be performed after neoadjuvant chemotherapy to assess oligometastatic colorectal cancer (CRC) treated surgically with curative intent. We also investigated urine as an alternative analyte for ctDNA MRD detection in this nongenitourinary setting. PATIENTS AND METHODS: We applied AVENIO targeted next-generation sequencing to plasma, tumor, and urine samples acquired on the day of curative-intent surgery from 24 prospectively enrolled patients with oligometastatic CRC. Age-related clonal hematopoiesis was accounted for by removing variants also present in white blood cells. Plasma and urine ctDNA MRD were correlated with tumor cells detected in the surgical specimen, and adjuvant treatment strategies were proposed based on ctDNA-inferred tumor mutational burden (iTMB) and targetable alterations. RESULTS: Seventy-one percent of patients were treated with neoadjuvant chemotherapy. Tumor-naive plasma ctDNA analysis detected MRD at a median level of 0.62% with 95% sensitivity and 100% specificity, and 94% and 77% sensitivity when only considering patients treated with neoadjuvant chemotherapy and putative driver mutations, respectively. In urine, ctDNA MRD detection specificity remained high at 100%, but sensitivity decreased to 64% with median levels being 11-fold lower than in plasma (P < .0001). Personalized ctDNA MRD oncogenomic analysis revealed 81% of patients might have been candidates for adjuvant immunotherapy based on high iTMB or targeted therapy based on actionable PIK3CA mutations. CONCLUSION: Tumor-naive plasma ctDNA analysis can sensitively and specifically detect MRD in patients with oligometastatic CRC after neoadjuvant chemotherapy. Urine-based ctDNA MRD detection is also feasible; however, it is less sensitive than plasma because of significantly lower levels. Oligometastatic patients with detectable MRD may benefit from additional personalized treatment based on ctDNA-derived oncogenomic profiling.

[Potential geographical distribution and changes of Artemisia ordosica in China under future climate change]. / 未来气候变化下黑沙蒿在中国的潜在地理分布及变迁.

Artemisia ordosica is a forerunner species of wind-break and sand-fixation in desert steppe in China, which plays an important role in ecosystem restoration and reconstruction. How-ever, it could influence human health. Based on 89 valid data of current distribution of A. ordosica in China and 19 typical climatic factors, the MaxEnt model was used to simulate the potential distribution of A. ordosica in China under current and two scenarios (RCP 4.5 and RCP 8.5; 2050s and 2070s). The SDM toolbox of ArcGIS software was used to analyze the potential distribution range of A. ordosica and its changes in China. The importance of key climatic factors was evaluated by comprehensive contribution rate, Jackknife method, and response curve of environmental variables. The accuracy of model was tested and evaluated by area under the curve (AUC) of the test subject working characteristic (ROC). The results showed that the MaxEnt model worked well (AUC=0.980). which predicted that A. ordosica was mainly concentrated in and around Mu Us Sandy Land, consistent with the current actual distribution range. The distribution area of A. ordosica of potential high fitness under the future two scenarios decreased by 5.2%-26.8%, which was negatively affected by future climate change. Seasonal variation of temperature, mean precipitation in the coldest season, and mean annual temperature had the greatest impact. The core area of future potential distribution of A. ordosica in China was located in Mu Us Sandy Land, with a tendency for spreading to northeast (Jilin, Heilongjiang, Liaoning and some parts of Hebei).

Cell-free DNA alterations in the AR enhancer and locus predict resistance to AR-directed therapy in patients with metastatic prostate cancer.

PURPOSE: Cell-free DNA (cfDNA) and circulating tumor cell (CTC) based liquid biopsies have emerged as potential tools to predict responses to androgen receptor (AR)-directed therapy in metastatic prostate cancer. However, due to complex mechanisms and incomplete understanding of genomic events involved in metastatic prostate cancer resistance, current assays (e.g. CTC AR-V7) demonstrate low sensitivity and remain underutilized. The recent discovery of AR enhancer amplification in >80% of metastatic patients and its association with disease resistance presents an opportunity to improve upon current assays. We hypothesized that tracking AR/enhancer genomic alterations in plasma cfDNA would detect resistance with high sensitivity and specificity. METHODS: We developed a targeted sequencing and analysis method as part of a new assay called Enhancer and neighboring loci of Androgen Receptor Sequencing (EnhanceAR-Seq). We applied EnhanceAR-Seq to plasma collected from 40 patients with metastatic prostate cancer treated with AR-directed therapy to monitor AR/enhancer genomic alterations and correlate these events with therapy resistance, progression-free survival (PFS) and overall survival (OS). RESULTS: EnhanceAR-Seq identified genomic alterations in the AR/enhancer locus in 45% of cases, including a 40% rate of AR enhancer amplification. Patients with AR/enhancer alterations had significantly worse PFS and OS than those without (6-month PFS: 30% vs. 71%, P=0.0002; 6-month OS: 59% vs. 100%, P=0.0015). AR/enhancer alterations in plasma cfDNA detected 18 of 23 resistant cases (78%) and outperformed the CTC AR-V7 assay which was also run on a subset of patients. CONCLUSION: cfDNA-based AR locus alterations, including of the enhancer, are strongly associated with resistance to AR-directed therapy and significantly worse survival. cfDNA analysis using EnhanceAR-Seq may enable more precise risk stratification and personalized therapeutic approaches for metastatic prostate cancer.